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Abstract 

Pulsar wind nebulae are now well established as important probes both of neutron 
stars' relativistic winds and of the surrounding interstellar medium. Amongst this 
diverse group of objects, pulsar bow shocks have long been regarded as an oddity, 
only seen around a handful of rapidly moving neutron stars. However, recent efforts 
at optical, radio and X-ray wavelengths have identified many new pulsar bow shocks, 
and these results have consequently motivated renewed theoretical efforts to model 
these systems. Here I review the new results and ideas which have emerged on these 
spectacular systems, and explain how bow shocks and "Crab-like" nebulae now form 
a consistent picture within our understanding of pulsar winds. 
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1 Introduction 



Pulsars have high space velocities, typical ly in the range 150-1500 km s ^ 
jArzoumanian. Cher'noff. and CorHe.j .m: Since such speeds are well in ex- 
cess of typical sound speeds for the cold (~ 1 km s~^) or warm components 
(~ 10 km s^^) of the interstellar medium (ISM), we expect that many pul- 
sars are moving supersonically, and that they should thus drive bow shocks 
through ambient gas. 



Such phenomena are, of course, not unique: astrophysical bow shocks are seen 

around proto-st ars, clumps of supernova ejecta, runaway OB stars, and in the 

solar wind (e.g., Kaper et al. , 1997 ; van Buren and McCrav , 19881 : Sutherland and Dopita 
199,4 IZankl . [l99£ l). However, several factors make pulsar bow shocks particu- 
larly appealing to study. First, for pulsars the mechanical luminosity of the 
wind driving the bow shock can be reasonably estimated through the pulsar's 
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"spin-down luminosity", given hj E = 4tt'^IP/P^ (where / = 10^^ g cm^ is 
the neutron star's assumed moment of inertia, P is the star's spin period, and 
P is the period derivative). Spin-down luminosities for observed pulsars are 
typically in the range 10^^ — 10^^ ergs s~^. Second, for many pulsars, precision 
timing and interferometric measurements allow accurate determinations of the 
star's position and proper motion. Finally, compared to many other Galactic 
populations, the distances to pulsa rs are reasonably well-e stablished — from 



parallax for about 20 sources (e.g.. IChatterjee et al 



I. .20041) . and from disper- 
Cordes and T;azinl. Eo02V 



sion of the pulsed radio signals for over 1500 others 
All these factors mean that many of the unknowns generally associated with 
bow shocks are not an issue for pulsars; in principle, the only parameters to 
be determined are the density and structure of the ISM, plus the inclination 
angle of the system. 

Furthermore, most pulsars are 10^ — 10^ years old, their high space velocities 
having carried them far from their birth places. This means that pulsars are 
more or less randomly located within the Galactic disk, and that their bow 
shocks can thus act as unbiased probes of the ISM. Bow shocks also potentially 
provide an opportunity to study the winds of "typical" pulsars, rather than 
extreme cases of youth and energy as typified by the Crab Nebula. 



Despite their potential as useful probes of pulsars and their surroundings, very 
few pulsar bow shocks were known until recently. Because of several recent 
discoveries and accompanying theoretical modelling, there has been renewed 
interest in these systems in the last few years. Here I summarise recent work, 
and the resultin g impro v ement s to our understanding. For e arlier reviews on 



bow shocks see CordesI fll99fit ) and Chatter jee and CordesI ()2002t ): for more 



general d i scussi on of pulsar wind nebul ae, the reader is referred to reviews by 
Gaensleii ()2004t ) and lKaspi etHI (|200.'5h . 



2 Theoretical Expectations 



Several authors have co nsidered the theory of bo^y shocks in d etail, both under 

996[lWilkinl.E99fi[l200' 



general circumstance s (iChen et al. 



a nd specifically 

I992I : Bucciantinil . l2002al lbl:l Ivan der Swaluw et al.l 



in the case of pulsars (jAldcroft et al 
2OO3I ). Some pertinent points from these studies follows. 



Scaling Parameter: One parameter sets the overall scale of a bow shock — 
the "stand-off distance" , r^. This is the separation between the pulsar position 
and the apex of the bow shock, and is set by pressure balance. 
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Fig. 1. Hydrodynamic simulation of a pulsar bow shock, adapted from Bucciantinil 
( 2nn2al ). The greyscale indicates density; the pulsar is moving from right to left with 
a Mach number M ~ 13. Various regions and interfaces discussed in the text are 
indicated. 

where p is the ambient mass density and V is the pulsar's space velocity 
relative to its surroundings. 



Shape: IWilkin has provided an elegant analytic solution to the shape 

of a bow shock, 



r{9) = CSC 0^/3(1 - ^cot 



(2) 



where r{6) is the radius of the bow shock at a polar angle 6. This solution 
is for an idealised thin-layer shock, but is also expected to be a reasonable 
approximation to pulsar bow shoc ks in regions near the apex (jBucciantini . 
2002a : van der Swaluw et al. . 2003[ ). 



Overall Structure: Because of the long cooling time scales, bow shocks are 
not thin shells, but show a characteristic double-shock structure, as shown 
in Figure 1. The morphology consists of a forward shock which heats the 
ISM, a termination shock which decelerates the pulsar wind, and a contact 
discontinuity dividing shocked ISM from shocked pulsar wind material. 

Termination Shock: The pulsar wind experiences a significant difference 
in pressure between regions ahead of and behind the pulsar's motion. The 
termination shock is therefore not of uniform radius around the pulsar, but 
rather has a bullet-shaped morphology, the head of the bullet aligned with 
the pulsar's direction of motion as can be seen in Figure 1. For low Mach 
numbers, ~ 1, the ratio of termination shock radii between the directions 
immed iately behind and directly ahead of the pulsar i s roughly proportional 
to M. ( Bucciantini . 2002a; van der Swaluw et al. . 20031 ). However, in the limit 
of high Mach number, this ratio tends to an asymptotic limit of about 5 
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f Gaensler et al 



20041 ) ■ 



3 Observations: Forward Shock 



The forward shock should produce observable Ha emission, resulting from 
collisional excitation and charge exchange of neutrals in ambient gas. Indeed 
some beautiful bow shocks have been observed in the Ha line, several of which 



have been discovered iust in the last few yea rs ()van Kerkwiik and Kulkarni , 
200ll : I.Tones et all I2OO2I : lOaensler et Zl . l2002h . Detection of such a bow shock 



allows one to estimate measure r^, although two cautions must be issued. First, 
formally corresponds to the distan ce from the pulsar t o the termination 



shock, not to the forward shock (e.g., Bucciantinl 2002a). So the observed 



pulsar/shock separation needs to be scaled by some factor, which simulations 
sugge st has a value ~ 0.4 — 0.6 ()Bucciantini . 2002a; van der Swaluw et al.L 
I2OO3I ). Second, the measured separation needs to be corrected for the (usually 
unknown) inclination. This correction is not a simple trigonometric factor, 
because projection of the three-dimens ional bow-shock surface onto the sky 
can be larger than the true separation ( Gaensler et al. . 20021) . Correction for 
inclination is non-trivial, and quantities inferred from the observed bow-shock 
separation should thus be considered approximate. 



With these caveats in mind, pressure balance can be applied to a system in 
which E and V are known to yield p via Equation (1). This approach suggests 
ambient densities ~ 0.1 cm~^ , as expected for warm neutral gas in the ISM 
fIChatteriee and Cordesl . I2OO2I : lOaensler etaP . l2002> ). . n cases where V is not 
known, one can write pV^ = 'yAi^V, where 7 = 5/3 is the adiabatic coefficient 
of the ISM and V is the ISM pressure. Adoption of a typical value for V allows 
one to directly estimate M.. 



In the case of an isotropic pulsar wind propagating through a homogeneous 
ISM, Equation (2) should describe the Ha brightness profile well. For PSR J0437- 
4715 this indeed seei ns to be the case, a llowing one to solve for the inclination 
angle of the system ()Mann et al. . 19991 ). However, in several other cases there 
are significant deviations from this idealised solution, in the form of abrupt 
kinks and bulges in the bow-shock profile, or a rotational offset between the 
symmetry axis o f the bow shock a nd the velocity vec tor of the pulsar ( e.g., 



PSR B0740-28, I.Tones et al.l I2OO2I : PSR J2124-3358, lOaensler et Zl I2OO2. ) 



These shapes imply some combination of a nisotropies in t he pulsar wind (as 
are known to exist around young pulsars; see lGaenslei]l2004[ ). density variations 
in the ISM, or significant motion of ambient gas with respect to the local rest- 
frame. Multi-epoch observations, showing changes in the bow-shock structure 
with time, are beginning to explore these possibilities ( Chatteriee and Cordesl . 
20041 ). 
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4 Observations: Termination Shock 



4.1 The Mouse 



Just as is the case for younger systems hke the Crab Nebula, we expect that 
the particles in the pulsar wind will be accelerated at the termination shock, 
generating radio and X-ray synchrotron emission. Indeed cometary radio and 
X-ray emission aligne d with the direction o f moti on is seen around several 
high velocity pulsars ( Chatteriee and Corded . 2002[ ). 



Figure 2 shows the results of a recent detailed study of "the Mouse", an 
elong ated radio and X-ray nebula coincident wi th the energetic pulsar J 1747- 
2958 (|Camilo et all 120021 : iGaensler et al.l . 12004^ . The Mouse is very luminous. 



and thus represents an ideal opportunity to test the theories developed for 
bow-shock termination shocks. 

If the X-ray and radio emission from the Mouse comes from the termination 
shock region, we expect it to be sharply bounded at its apex by the contact 
discontinuity. The extent of emission thus allows us to estimate the stand-off 
distance for this system — for a distance of 5 kpc we find ~ 0.02 pc. For 
ISM pressures V/k = 2400Po K cm~'^ (with typical observed values in the 

— 1/2 

range 0.5 J: Pq ^ 5), this implies a high Mach number, A4 ~ SOPg > cind 
thus a likely space velocity V ~ GOOPg km s^^ if propagating through the 
warm component of the ISM. 

The X-ray and radio morphologies shown in Figure 2 both indicate the pres- 
ence of an elongated bright "tongue" of emission behind the pulsar, about 
0.25 pc in length. An interpretation of this feature is suggested by compar- 
ison with the hydrodynamic simulation shown in Figure 1, where it can be 
seen that the tongue has a similar appearance to the surface of the bullet- 
shaped termination shock. If we make this identification, then the tongue in 
the Mouse is directly analogous to the inner toroidal ring seen for the Crab 



Nebula and other related systems (e.g.. iNg and Romanil . |2004( ). demarcating 
the point where the pulsar wind reacts to its surroundings and where wind 
particles are accelerated up to synchrotron-emitting energies. 

It is important to note that hydrodynamic simulations such as that shown in 
Figure 1 do not explicitly include the effects of magnetic fields or relativistic 
flows, so an exact correspondence between the data and the model is not ex- 
pected. In particular, two important discrepancies can be identified: the ratio 
of backward and forward termination shocks is observed to be ~ 10, rather 
than the value of 5 predicted for ^ 1; and the sheath does not show 
any limb-brightening as would be expected given its geometry. The first of 



5 



o 
o 
o 

M 



U 
H 

a 



48 
54 
58 00 
06 
12 - 
-29°58'l8" - 



(a) Chandra (0.5-8,0 keV) 




tongue 



tail 



o 
o 
o 



u 

Q 



48" : (b) VLA (4.8 GHz) 
54" - 

58 00" - 
06" - 

12" - narrow 

-29°58'i8" - broad 




o 
o 



o 

CO 



o 

CD 



O 



o 

C\2 



I 



17''47™17" 



IB"- 



15^ 14^ 
RA (J2000) 



13" 



Fig. 2. X-ray (upper) and radio (lower) images of the bow shock G359 . 23-0.82 



(a.k.a. "the Mouse") powered by PSR J1747-2958 (|Gaensler et al.l . l2nn4h . In the 



X-ray image, the "tongue" and "tail" regions are indicated. In the radio image, the 
broad and narrow components of the tail are shown. 

these points can potentially be explained by the effect of ions in the wind, 
whose gy rations in the magnetic fi eld spread the electron shock over a large 
distance ( Gallant and Arona . 1994 ). This effect should operate in the back- 
ward flow from the pulsar, but not in the forward direction due to the high 
degree of pressure confinement. The lack of limb brightening can be accounted 
for by Doppler boosting in the post-shock flo w', as is consistent l y obs erved in 
the winds of younger, lower velocity pulsars ()Ng and Romanil . 12004^ . While 
conflrmation of these effects will require a fully relativistic magnetohydrody- 
namic treatment, we emphasise that both the effects claimed here have been 
seen in the wind nebulae around "Crab-like" pulsars. The need to invoke sim- 
ilar phenomena here suggests a high degree of commonality between the wind 
properties of the youngest pulsars and more typical members of the popula- 
tion. 



West of the tongue, a fading tail of emission is seen, which must come from 
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material significantly downstream of the termination shock. This emission 
is analogous to the main body of the nebula in Crab-like systems. The tail 
of emission behind the Mouse can be delineated into two components, as 
indicated in Figure 2. The first component is bright, reasonably narrow, and 
is seen in both radio and X-rays; the second component is broader, fainter, 
and is seen in radio only. These structures can be qualitatively understood by 
again considering the hydrodynamic simulation of Figure 1, in which it can 
be seen that there are two extremes in the flow. Ahead of the pulsar, the wind 
shocks at a small distance from the pulsar at the apex, then turns around and 
flows around the edges of the shocked region. In contrast, behind the pulsar the 
distance to the shock is considerably larger; the post-shock flow then remains 
in a narrow collimated region lying along the symmetry axis. If equipartition 
holds, so that E/Airr'^c oc B^, where B is the post-shock nebular magnetic 
field, then the shocked fiow ahead of the pulsar should have a magnetic field 
~ 10 times larger than that behind the pulsar. We thus expect a broad fiow 
in which the X-ray emission suffers severe synchrotron losses, enveloping a 
narrow flow in which the X-rays are yet to cool. This is in accordance with 
observations, suggesting that in contrast to the Crab Nebula, the pulsar's 
motion structures the post-shock emitting regions into two distinct zones. 



4-2 X-ray Trails Behind Other Pulsars 



Now that we have a picture of what might be occurring for the Mouse, it 
is insightful to consider other bow shocks which produce radio and X-ray 
emission, to see if similar features are observed. 



Some p ulsars such as PSRs B1757-24 and B1957+20 (jKaspi et al.L 120011 
Istapper s et all l2003t) , show short (~ 1—0.5 pc) X-ray trails extending oppo- 
site the pulsar's direction of motion. One interpretation is that these features 
are synchrotron "wakes", i.e., particles left behind by the pulsar as it moves 
through the ISM. However, the time taken for a pulsar to traverse the extent 
of one of these trails is ~ 1000 yr, implying abnormally low nebular magnetic 
fleld strengths if synchrotron cooling at X-ray energies is to be avoided during 
this period. It has therefore been suggested that a rapid backflow or nozzle op- 
erates downstream, which rapi dly advects emitting particles and produces th e 



erates aownstream, wnicn rapi dly advects emittmg particles and produces tm 
trail structures observed (e.g., Wang and Gotthelf . Il998l : Kaspi et al. . 2001). 



However, the Mouse suggests an alternative interpretation. We expect two 
components to the downstream synchrotron emission from a pulsar bow shock. 
Close to the pulsar, we should see a "tongue" of emission around the termina- 
tion shock. We expect this region to show a constant X-ray spectrum across 
its extent, to have a comparable extent in radio and X-rays, and to abruptly 
terminate at a distance ~ lOr^ behind the pulsar. Further from the pulsar. 
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Fig. 3. X-ray and radio observations of tlie X-ray bow shock around 
CXOU J061705.3-F222127 in the SNR IC 443. The image shows Chandra data 
smoothed with a 2" gaussian, while contours denote 8.4-GHz VLA data at a reso- 
lution of 7". Proposed features are indicated. 

the "tail" component representing post-shock material should show an X-ray 
spectrum which softens with distance from the pulsar, should have a larger 
extent in radio than in X-rays, and should gradually fade into the background 
over an extent ^ lOr^- 



Armed with the realisation that bow shocks should show both these compo- 
nents, it seems that the short trails seen around PSRs B 1 957-1-20 and B 1757-24 
corr espond well to the to ngue component of the Mouse. iGaensler et al.l (|2004 ) 
and Gvaramadzd ( 2004 ) have thus both argued that the name "trail" is a 
misnomer — these features are simply the surfaces of the termination shock 
around these pulsars, the longer tail of emission downstream being too faint 
to see. 



On the other hand, a closer examination of other pulsar bow shocks suggests 
that they too potentially show multiple components in X-rays, corresponding 
both to the termination shock and to the tail of emission behind it. Possible 
exam ples of such systern s are PSR B1853+01 in the supern ova remnant (SNR ) 
W44 (|Petre et al.Ll2002h . PSR B1951+32 in SNR CTB 80 (iMoon et al.Ll2004l) 
and t he neutron star CXOU J061705.3+222127 in SNR IC 443 (|01bert et all 



20011 ). The latter source is shown in Figure 3, where the identification of these 
possible structures is indicated. It is interesting to note that this source is 
moving through the shocked gas in the interior of an evolved SNR. In this 



case, we expect a Mach number ~ 3 ( van der Swaluw et al. . 2003[ l. much 



lower than for a pulsar moving through the ISM. For this low Mach number, 
we expect the "tongue" region to only be about half as long (in units of 
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Ty^) than for a high Mach number system hke the Mouse. Indeed the data in 
Figure 3 suggest an extent for this feature ~ 5r^, compared to ~ lOr^ for 
the Mouse. Deeper Chandra observations of both CXOU J061705.3+222127 
and PSR B1853+01 are scheduled to take place in 2005, which should allow 
a better investigation of these possibilities. 



5 Putting It All Together 



In §2, it was noted that a double-shock structure is expected around a su- 
personic pulsar. While we clearly see either forward or termination shocks 
around various neutron stars, for just one system, PSR B 1957-1-20, have both 
the termination sho ck (in X-rays) and for ward shock (in Ha) been identified 
in the same source ( Stappers et al. . 2003| ). A collection of such sources can 



act as a detailed probe of the hydrodynamics of the bow-shock interaction, so 
expansion of this sample is obviously highly desirable. 

Most of the pulsar bow shocks seen in radio/X-rays have high values of E, 
meaning that they are reasonably rare and hence distant. The high extinction 
to these systems makes it difficult to identify any Ha emission around them. 
On the other hand. Ha bow shocks tend to be nearby, so are worthwhile targets 
for X-ray imaging, even if E is low. But of the five other kno wn Ha bow shocks 
besides that around PSR B1957+20, two (PSR .10437-4715. fZavhn et allliooi 
RX J1856.5-3754 M show no extended X-ray emission, one (PSR B2224+65; 
Wong et al 1l2003t) shows an X-ray filament completely misaligned with the pul- 



sar's direction of motion, and the remaining two (PSRs B0740-28 and J2124- 
3358) are yet to be observed at high angular resolution by Chandra. Radio 
searches to some o f these sources have a l so not been successful at det ecting 
extended emission (|Gaensler et al.l . I2OOOI : IChatteriee and Cordesl . l2002> ). The 



lack of emission from the termination shock in systems for which the wind is 
known to be strongly confined is puzzling. 



6 Conclusions 



A consistent understanding is now beginning to emerge of the various multi- 
wavelength structures seen around high velocity pulsars. As expected, we see 

^ This source is one of several nearby isolated neutron stars seen in thermal X-rays, 
and does not pulse. We assume that this source is a typical pulsar but beaming away 
from us. To the best of our knowledge, the lack of extended X-ray emission has not 
been stated in the literature, but is immediately apparent from inspection of deep 
archival Chandra observations of this field. 
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optical emission from the shocked ISM in several bow shocks, from which 
we can infer the density of the ambient medium or the Mach number of the 
system. The morphologies of these shocks also allow us to probe anisotropics 
in the ISM and in the pulsar wind. 

The X-ray and radio emission produced at the termination shock in these 
systems permit a study of the relativistic wind which drives the bow shock. 
Just as in the Crab Nebula, high-resolution X-ray imaging allows us to identify 
both the surface of the termination shock and emission from the post-shock 
flow. As in the Crab, the effects of an ion-loaded wind and of Doppler beaming 
are required to explain the observations, suggesting that the wind composition 
and behaviour is similar for bow shocks and for Crab-like systems. However, a 
distinct difference from the Crab is the presence of two flow zones downstream, 
corresponding to highly- and weakly-magnetized material, emerging from the 
front and back of the termination shock, respectively. 

In the near future, we expect that deeper X-ray and opti cal images will 



increase the sample and show new types of structure (e.g., ICaraveo et al. 



2003| ). and that relativistic magnetohydrodynamical simulations can begin to 



provide a solid theoretical underpinning for the observed morphologies. The 
time domain is now also beginning to play a role, in that we can watch bow 
shocks evolve as they trace the cascade of density inhomogeneities in the ISM 
f Chatteriee and Cordesl . 2004 ). Clearly bow shocks are now beginning to ful- 



fil their long-standing promise to act as unique probes of pulsars and their 
surroundings. 
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